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ABSTRACT 

Gamma-ray burst statistics are best explained by a source population at cosmological distances, while spec- 
troscopy and intensity histories of some individual bursts imply an origin on Galactic neutron stars. To 
resolve this inconsistency I suggest the presence of two populations, one at cosmological distances and the 
other Galactic. I build on ideas of Shemi & Piran (1990) and of Rees & MesSz&ros (1992) involving the inter- 
action of fireball debris with surrounding clouds to explain the observed intensity histories in bursts at cosmo- 
logical distances. The distances to the Galactic population are undetermined because they are too few to affect 
the statistics of intensity and direction; I explain them as resulting from magnetic reconnection in neutron star 
magnetospheres. An appendix describes the late evolution of the debris as a relativistic blast wave. 

Subject headings: gamma rays; bursts — MHD — stars: neutron 


1. INTRODUCTION 

Attempts to explain all the observed gamma-ray bursts 
(GRBs) with a single population of sources have become pro- 
gressively more difficult. On one hand, their distribution on the 
sky has been observed, with steadily improving precision 
(Atteia et al. 1987; Megan et al. 1992), to be isotropic, an 
observation which is naturally explained (Usov & Chibisov 
1975; Goodman 1986; Paczynski 1986; Dar et al. 1992; 
Dermer 1992; Fenimore et al. 1992; Mao & Paczynski 1992a; 
Meszaros & Rees 1992; Paczynski 1992a; Piran 1992) if they 
are at cosmological distances. On the other hand, a number of 
GRBs have been reported to show (Higdon & Lingenfelter 
1990) spectral features at a few tens of keV and at about 400 
keV, which are readily interpretable as cyclotron lines and the 
two-photon positron annihilation line from the surface of mag- 
netized neutron stars at Galactic distances, but which are inex- 
plicable at cosmological distances. If their validity is accepted, 
the data appear irreconcilable. 

The problem is further complicated by the fact that straight- 
forward models of radiation transport in GRBs at cosmo- 
logical distances (Goodman 1986; Paczynski 1986) predict 
very brief bursts of radiation with thermalized spectra, in con- 
tradiction to observation, while attempts (Brainerd 1992; 
Harding & Leventhal 1992; Katz 1992; Li & Dermer 1992; Li 
& Liang 1992; Mao & Paczynski 1992b; Wasserman 1992) to 
explain the spatial anisotropy and log N versus log 5 or 
V/V mtkX distributions of GRBs in Galactic models require the 
assumption of a spherically symmetric halo of ^100 kpc 
radius. Finally, the soft gamma repeaters (SGRs) introduce 
additional confusion. The fact of their repetition and the iden- 
tification of one of them (1979 March 5; Cline 1980) with a 
supernova remnant in the LMC point strongly to a Galactic 
population, while the presence of spectral features and an 8 
second periodicity (1979 March 5) indicate origin on a mag- 
netic neutron star. However, it is unclear that SGRs should be 
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considered GRBs at all because their properties, jnduding 
their spectra, are very different, and arguments made for SGRs 

may be irrelevant to the problems of GRBs. 

As a first step toward resolving the apparent inconsistencies, 

I consider the obvious possibility that there are two distinct 
populations of GRBs. A cosmological population C accounts 
for most GRBs, and explains the statistics of isotropy and the 
log N versus log S or V/V ml distributions. A population G 
includes those GRBs with spectral lines (a minority), and orig- 
inate on neutron stars at Galactic distances. An individua 
GRB cannot be assigned to a population unless it shows spec- 
tral lines, but the majority (probably overwhelming) of i 
without spectral lines must be members of population C. lne 
SGR may be members of population G. I do not consider more 
exotic possibilities, such as GRBs arising within the Oort cloud 
(Ruderman 1975), because they ail seem fa r ' fet ^J 
though I myself have discussed one of them (Katz 1993). Th 
suggestion of a cosmological population (to explain the 
statistics) and a Galactic population (to explain the spectral 
lines) is not related to the suggestion of Lingenfelter & Higdon 
(1992) (disputed by Paczynski 1992b) that the statistics may be 
explained by two Galactic populations. 

In the absence of information about the intensity statistics 
and angular distribution of population G alone, it is not pos- 
sible to discriminate between a disk and a halo origin. It will be 
important to obtain this information, which may be reducible 

from archival data. r __ _ 

In § 2, 1 discuss a possible mechanism for GRBs at cosmo- 
logical distances, building on recent suggestions by others. In 
6 3 I more briefly discuss magnetic reconnection models oi 
GRBs at Galactic distances. The 1979 March 5 event in the 
LMC poses an acute problem of gamma-gamma pair pro- 
duction (Carrigan & Katz 1992) which must be faced, whether 
or not GRBs of population G have comparable distances and 
luminosities, and even if there is no population G. § 4 contains 
a summary discussion. Unfortunately, unambiguous observa- 
tional tests of the ideas discussed here will not be easy. 

2. POPULATION C: GAMMA-RAY BURSTERS AT 
COSMOLOGICAL DISTANCES 

The well-known failure of straightforward fireball models to 
explain the spectral and temporal properties of GRB led Shram 
& Piran (1990) to consider neutrino-produced fireballs loaded 
with small (but not zero) amounts of ordinary matter; they 
found that the fireball could couple nearly all of its energy to 
the matter and (with the right values of the parameters) could 
accelerate it to relativistic velocity. Rees & Meszaros (199ZJ 
and Meszaros & Rees (1993) then pointed out that the inter- 
action of this relativistic debris with surrounding matter such 
as that suggested by Narayan, Paczynski, & Piran (1992), 
might be characterized by times consistent with the range ot 
GRB rise times and durations of(10' 2 3 -10 3 s). 

I build on these ideas. The environments of GRBs at cosmo- 
logical distances are open to much speculation (for example, 
are they low-density Galactic halos or dense nuclei of gal- 
axies?), but the strong dumpiness of interstellar matter is a 
consequence of immutable atomic physics (cooling rates), and 
isolated discrete clouds are likely under a very wide range of 
conditions. The rarity of GRBs makes it possible to assume 
favorable conditions, if these lie in the range of plausibility, 
there is no great difficulty if a considerably larger number ot 
fireballs occurring in less favorable rircumstances do not 
produce observable GRBs. 

GRBs at cosmological distances require the radiation of 
~ 10 51 ergs in observable gamma-rays. The complex chain ot 
processes which lead to gamma-ray emission must be moder- 
ately efficient when the parameters do have favorable values, 
because the energy radiated as neutrinos by neutron star col- 
lapse, formation, or coalescence is unlikely to exceed 0.3 M© 


C 2 « 5 X 10” ergs, and may be considerably less. Thus, while 
we are entitled to assume favorable circumstances to explain 
the rare observable GRBs, when these circumstances occur the 
resulting processes must be reasonably efficient. The fraction of 
neutrino energy converted to an electromagnetic fireball is 
small. Efficient conversion requires neutrino-neutrino colli- 
sions at angles in excess of 90°, but the neutrinos are generally 
expanding outward in a neutrino fireball, with velocity vectors 
which are tending toward radial outflow. Optimal head-on 
collisions are particularly rare. The conversion of electromag- 
netic energy to particle kinetic energy also has an efficiency less 
than 1. The final conversion to observable gamma-rays is the 
hardest part of the problem; it is easy to see how this could fail 

entirely. , 

Relativistic invariants alone limit the amount of kinetic 
energy available for radiation by fireball debris. If a relativistic 
debris cloud with speed /J f c, Lorentz factor y F , and proper 
mass per unit area a sweeps up a proper mass per unit area a<r, 
then the efficiency of radiation of the debris kinetic energy can 
be as large as 


a + VrO-W 

Values of e > i are obtained for a > l/[v^<l + Ml « l/(2y F ). 
Collisions with a very broad range of clouds of circum-fireball 
matter are consistent with efficient conversion of kinetic energy 
to radiation. This is fortunate, because efficient production of 
GRBs requires that common circum-fireball environments 
produce observable GRBs in most directions; it is not possible 
to insist on fortuitous geometries or on special values of the 
parameters. 

Most of the kinetic energy will become available when the 
debris has swept up only a very little matter. If the energy of 
the explosion is Y then the proper mass of debris is Y/y F c 2 , 
and half the kinetic energy will become available when the 
swept-up proper mass is Y/2y F c 2 . In a uniform medium of 
density p this will occur at an interaction radius 



~ 2 x 10 15 cm , 


( 2 ) 


where the numerical estimate assumed Y = 10” ergs, 
p = 10' 24 gem' 3 and y F = 10*. 

The hardest part of the problem is turning the kinetic energy 
of the relativistic debris into the observed gamma-rays. The 
collision length of relativistic protons in ordinary matter is 
about 50 g cm' 2 , or about 100 Mpc at typical interstellar 
densities. Clearly, some collective process is necessary, and it 
must couple the proton and ion energy into that of electrons, 
which radiate more readily. Even relativistic electrons do not 
radiate rapidly under interstellar conditions; the radiation 
length of a 10 13 eV election (corresponding to equipartition 
with a y F ~ 10 4 proton^) for Compton scattering on a 3 K 
blackbody radiation field is ~10 23 cm, excessive by many 
orders of magnitude. 

In order to obtain short pulses of radiation at distances of 
order those given by equation (2) it is necessary that a coherent 
relativistically expanding front of radiating particles be 
directed nearly toward the observer. It is not sufficient that 
individual particles be observed only when directed toward the 
observer, a condition met by most relativitlstic radiation pro- 
cesses. Therefore, ambient magnetic fields must not deflect par- 
ticles significantly from their initial spherical expansion. This 
condition will be satisfied if the magnetic energy E Mt in the 
interaction sphere of radius r, is very much less than the debris 
energy so the debris can sweep away the ambient magnetic 
field without significant deflection. If equipartition is assumed 
between the ambient magnetic field B|s M and an ambient turb- 
ulent velocity field t) ISM , then E M JY ~ v ls Jy F c < l/y f , so 


that the ambient field may safely be ignored. 

2.1. Shock Structure 

When the debris shell collides with a cloud of ambient 
matter the resulting flow may be complex. If the shell and the 
cloud each initially had uniform density and velocity and negli- 
gible (on a relativistic scale) temperature, the geometry is slab- 
symmetric, and all bulk velocities are normal to the planes of 
symmetry, then the resulting shock structure is shown in 
Figure 1. There are two shocks SI and S2 and, in general, a 
contact discontinuity CD separating shocked fireball debris 
from shocked cloud. The equations relating the conditions in 
the four regions are cumbersome, except in the special sym- 
metric case in which debris and cloud initially had the same 
composition and proper density. In this case, which I assume, 
there is no contact discontinuity and conditions in regions 2 
and 3 are identical, as are those in regions 1 and 4. 


Fig. 1. — Flow geometry in frame of contact discontinuity CD. SI and S2 
arc shocks and numbers denote regions of fluid. 


The relativistic shock conditions (Landau & Lifshitz 1959) 
may be used to determine physical conditions. The thermody- 
namic variables are the proper internal energy density e, the 
proper pressure p, and the proper enthalpy density w = e + p. 
In the unshocked cloud e x = where is the proper 

atomic number density and m a is the proper mass per atom, 
and Pi = 0. In the shocked cloud p 2 = e 2 /3> in the 
extreme-relativistic (ER) limit. I shall refer to the frame of the 
unshocked interstellar material as the local observer’s frame; 
transformation to our frame requires application of the cosmo- 
logical redshift. Then, to lowest nontrivial order in eje 2 < 1, 
the velocities of the fluids with respect to the shock front SI are 


\ — — 

e 2 


and 1 


V± _ Rpa ~ PiXea + Pi) ! 1 
c |_(e 2 - e,Xei + P 2 )J 

Vi _ f (P 2 ~ PiX?i + Pz) ~] 12 as - f 1 + 2 — Y 
c L(e 2 - ei\e 2 + Pj) J 3 \ ej 


(3) 

(4) 


1 Note the assertion in the first edition of Landau and Lifshitz (1959) that in 
the ER limit v 2 -+ c/3 1 ' 2 is a typographical error; the correct limit v 2 -*■ c/3 is 
given in later editions. 

The velocity discontinuity v l2 between fluids 1 and 2, mea- 
sured in the frame of either, is obtained from the relativistic 
expression for the subtraction of velocities 

V 12 _ gi / c ~ V 2 ! c ^ j _ 2 £l • (5) 

c 1 - v x vjc 1 e 2 

this velocity is also the velocity v 2L of shocked fluid 2 in the 
local observer’s frame. The velocity v x is also the speed of the 
shock SI in that frame. 

Fluids 2 and 3 have the same velocity and, given our 
assumptions that n 1 — rt 4 and e ^ = e 4 , the same values of the 
thermodynamic variables. Then the velocity of fluid 3 in the 
frame of shock S2 is - v 2 , and the speed of fluid 4 in that same 
frame is -u t . The expressions, for combinations of relativistic 
velocities may then be used to obtain the following results in 
the local ofcferver’s frame: 



c 


( 6 ) 


(7) 


— « 1 — 4 — , 

.C e 2 


a* « i - 2( - 

c \e 2 


( 8 ) 


It is now possible to calculate e 2 from the debris Lorentz 
factor y, , defined in the local observer’s frame, using equation 
( 8 ): 


1 


It = 


[1 - (i >Jc) 2 Y 12 
e 2 & 2n l m a c 1 y F . 


e_i 

2e, ’ 


(9) 

( 10 ) 

The Lorentz factor y 2L of the shocked material in the local 
observer’s frame is obtained from i> 2L = i>, 2 , equations (5) and 
(10): 

TiL-ri' 1 . 

The detailed mechanics of the shock are obscure, but must 
be collisionless in order to form a shock at all. The shocked 
matter need not be in thermodynamic equilibrium. Heating ol 
the shocked material by plasma instabilities is the source of 
dissipation; the distribution functions of particle energies will 
not be (relativistic) Maxwellians, but are more likely to be 
power laws. The distribution of energy between electrons and 
ions is uncertain. In a highly relativistic shock, as we expect 
here electrons and ions are kinematically very similar 
(identical in the ER limit), so I will assume that the distnbution 
of particle energies is independent of species; roughly half the 
post-shock energy resides in electrons. Any neutral matter does 
not interact with the shock, so that the density n, refers only to 
the ionized component. At distances ~ r, (eq. [2]) the inter- 
stellar material will largely have been ionized by the flash ot 
radiation associated with the fireball or by collision with 

debris. , , 

Using the shock jump conditions for the proper enthalpy w, 
w, = n { m a c 2 , and the ER limit w 2 « 4e 2 /3 yields the proper 

atomic density 


* 2 (m‘c 2 ) 


l/i 


Define y 2 by the relation 


Vi = 


n 2 m a c 


( 12 ) 


(13) 


then the mean energy (in the frame of fluid 2) per particle is 
v, um c 2 , where jim. is the mean proper mass per particle. I- or 
pure ionized hydrogen M = 0.5, while for the usual cosmic 
abundances (fully ionized) n = 0.62. The mean Lorentz factor 
of an electron (in the frame of fluid 2) is 


7u = Vi 


F 1 . 


(14) 


The proper density n 2 (eq. [12]) may be rewritten, using 
equation (13), as 

n 2 &4n l y 2 • 

reproducing the result n 2 = 4n, for a strong but nonrclativistic 
shock (y 2 - 1) in a gas with an adiabatic exponent of 5/3. The 
Lorentz factor y 2 is found from its definition (eq. [13]), and 
equations (9) and (15): 

(16) 


..-ifer-tr 


In the local observer’s frame most of the particles are narTOwly 
collimated in the direction of the motion of fluid 2, and the 


typical Loren tz factor is larger than those given in equations 
(14) and (16) by a factor ~y 2L (eq. [11]). The angular width of 
collimation depends on the angular distribution of the particle 
momenta in the proper frame of fluid 2, which is unknown. If 
this is isotropic, then the locally observed angular width (for 
electrons as well as ions) 

0 o ~ylL l «y^ (it) 

Note that this is a much broader angular distribution than the 
locally observed radiation pattern from a single particle, whose 
Lorentz factor is ~y F (ions) or ~ y F mjm t (electrons). 

2.2. Time Dependence 

The geometry of radiation from an advancing spherical 
shock is shown in Figure 2. The distant (but cosmologically 
local) observer may first see a flash of radiation from the fire- 
ball, itself, whose arrival time is taken as t = 0. If the fireball is 
a consequence of the merger of binary neutron stars (Eichler et 
al. 1989), the initial pulse includes bursts of neutrino and gravi- 
tational radiation, as well as electromagnetic radiation. Their 
emission is essentially simultaneous, although their arrival may 
be affected by dispersion arising from plasma refraction, neu- 
trino rest mass (if any), etc. The initial electromagnetic flash is 
expected to have a thermal spectrum and to be extremely brief 
(<|10* 4 s) because of the small size of the fireball (Goodman 
1986; Carrigan & Katz 1992); if, as assumed here, the fireball 
energy is largely converted (Shemi & Piran 1990) to kinetic 
energy of debris this initial flash may be unobservably faint. 
However, if it is observed the time interval between it and the 
rest of the GRB/is an important constraint on the emission 
geometry. 


Fig. 2 . — Emission geometry 


Radiation emitted from a point (r, 0) on the expanding 
spherical shell arrives at the observer at a time 


t 


r(l — cos 0) 
c 



r(l — cos 0) r 

Z + -, 

c u 


(18) 


where v is the shell’s expansion velocity and the parameter 
(dimensionally but not physically a velocity) u = vc/{c — v). If 
the angular distribution of radiated intensity, measured in the 
local observer’s frame, is /( O ' ), where O' is the angle from the 
normal to the radiating surface, then the energy dE radiated by 
a patch of area dA is 


dE=f(9')dA . (19) 


Radiation directed toward the observer has 0' = 0. Using 
dA = 2nr 2 sin 0 dd and dt = r sin dO/c yields the observed 
power 


p (t, r) = — = 2 Kcrf(0 ) . 


( 20 ) 


The function J[6) is proportional to a convolution of the 
angular distribution of the radiating particles and their radi- 
ation pattern ; as previously discussed, the latter is expected to 
be narrower than the former. A plausible guess is then 

m *¥Tel' (21) 

where the angular width 0 O ~ )v^is essentially the same as 
that of the momentum distribution (eq. [17]). The pulse shape 
is then obtained from equations (20) and (21), using equation 


(is; io eliminate t/f 


r 


P(t, r) oc« 


r r 

I [(f - r/u)2c/r + 0g] 1 ~ u 


( 22 ) 


0, 


T 

t < - 
U 


where the approximation cos 0 ss 1 - 6 2 /2 has been used. The 
pulse shape P(t y r) has been plotted in Figure 3, where the _ 
dimensionless parameter t = Ict/rBl has been defined. 


i oSt- 


cf itfr h 




Fig. 3. — Radiation pulse emitted at radius r, plotted for — 20c/u. 
Dashed lines denote rise and fall for a cloud sharply circumscribed between 
- ( T i^o “ 2c/u) 112 and 6 2 = {t 2 0q - 2 cju) l/2 . Note that abrupt rise at t = 
is obtained for a uniformly filled medium, independent of the angular 
distribution of radiation, and does not require a cloud bounded in an angle. 
Units of the ordinate are arbitrary. 


The pulse form of equation (22) should be regarded only as 
an envelope, for the actual pulse shape will be modulated by 
the spatial distribution of matter which the debris shell sweeps 
up. One striking feature of equation (22) is its abrupt rise and 
pronounced time-skewness (Weisskopf et al. 1978), consistent 
with the observed rapid rises of some GRBs (Bhat et al. 1992; 
Fishman 1993). The characteristic width of this function is 


At 



r 

JfC ' 


(23) 


~ Use of r ~ r, (eq. [2]) and y F ~ 10 4 yields At - 10 s, the right 
order of magnitude for the duration of GRBs. Much longer or 
shorter At may be possible for plausibly different values of the 
parameters, particularly the cloud density, which is uncertain 
even to order of magnitude. 

The debris shell and shock propagate into a very heter- 
ogeneous medium. The effects of structure in 0 are shown by 
the dashed lines in Figure 3, which assume a cloud uniform in 
the range & x <£ 6 ^ 0 2 , with abrupt boundaries. A more realis- 
tic gradual density profile or shape would produce a gradual 
rise and decay ; the abrupt rise remains if (and only if) the cloud 
includes the line 0 = 0. Thus this abrupt rise is expected for 
some, but perhaps not all, GRBs, in accord with observations. 

A complete intensity profile of a GRB requires the integra- 
tion of equation (22) over r : 

m = ["V r)g(r)dr , (24) 

where the weighting function g(r) includes both the fact that 
the energy available for radiation falls off for r > r f and the 
effects of dumpiness of the ambient matter as a function of r. 
The observable region is a narrow half-cone of apical angle 
^d 0 : unless the scale of spatial structure is <r l 6 0 ~ 2 x 10 13 
cm, dumpiness is more likely to be apparent as a function of r 
than of 0, justifying the use here of equation (22) which ignored 
any dependence of density on 0. 

In the absence of spatial heterogeneity ^(r) may be taken to 
impose a cutoff at r % r 7 , so that 

rmin(r/,uf) 

0>(t) * I P(t, r)dr o c 


1 


wt) 


let + r[— (2c/u) + 0g] 


(25) 

The integral is elementary, but cumbersome, and of limited 
quantitative interest because of the artificiality of the assumed 
uniform density; the rise time is r,/u. Two possibilities should 






1 Por a shock SI propagating through a homogeneous 
F° r _ / ^ i ; _ 1/(2v f ) and u ft; 2y F c. Then, for 

medium v/c * 1 d _ 1Q i the n Se time is several 

r, «*» »? ‘rg, and 7 consistent with very 

en ? e nn e the other hand, the fireball debris propagates through 
Vnnm or vS low density intercloud medium with a speed 
& J. i (1 _ 1 so that u * 2yjc; impacts upon small d»s 

£.7 cloud, sca«»* -Mf,* c. mS 

SgSiET Thl ab7up. rise given by equation (22, and in 

^'integration over r introduces two broadenings one 
Ol?llT*c associated with the entire emission region of size 
° ( r' where the low intercloud density is appropriate, and 

Sous rarefactions and 

tivistic hydrodynamic calculations. 

2.3. Radiation 

Th* hardest oart of the problem is turning electron kinetic 
the observed radiation. Even though the typical 

rr on?' s "®s* 

sr 

therefore make the radical suggestion that the colhsionless 
shockproduces approximate equipartition between the mag- 
neS eneVgy density and the particle energy density The 
proper magnetic field B 2 and energy density in fluid 2 are then, 

using equation (10), 

Bl „ . .. „ „2 (26) 


^ = Ce 2 *2Cni7 f m 0 c 2 , 
87C 


where ( ^ i is a phenomenological parameter describing the 
approach to equipartition. The synchrotron energy loss tme 

for an electron with Lorentz factor given by equation ( 
assurnirurTio correlation between the direction of the elector,, 
momentum and that of the magnetic field, is then, taking pure 
hydrogen composition, 


tr2*C 


-1 


1 cm" 


If 


c 3 ' 2 1.1 x 10 7 s . 


(27) 


The radiating volume is moving toward the observer with a 
bulk Lorentz factor y 2L (eq. [11]), so that application of a 
Lorentz transform yields the local observers measured radi- 
ation time 

khl \ 

t ob, = V2L^~~p J 

lr2_^ r -i(Sh Y' y r 2 1.1 x 10 7 s . (28) 

= V 2 L \1 cm -3 / 




Jv Uh) 



- 


c 


For a plausible interstellar cloud density n, > 1 cm 3 and 
v, a 10 5 t h. may be a millisecond or less, as required by rapid 
time-structure in some GRBs (Bhat et al. 1992). This justifies 
the assumption, made implicitly in the discussion of GRB rise 
times and pulse lengths, that shock-accelerated electrons 
radiate instantaneously; properly, the pulse profiles predicted 
by equations (24) and (25) should be convolved with a 
broadening function which includes the radiation time, and 
which has a width t ob , in the local observer’s frame. 

The characteristic frequency of synchrotron radiation, mea- 
sured in the frame of fluid 2, is obtained from standard expres- 
sions using equautions (14) and (26). For pure hydrogen the 
result is 


\1 cm 


1/2 


yf 


3 x 10 u s' 1 , 


(29) 


while Lorentz transformation to the local observers frame, 
using equation (1 1), yields 


v 0 b» 


^ \1 cm’ 


1/2 


yj 3 X 10 u s 


(30) 


Mev photons may be observed for { - ± if n, ~ j cm 3 and 
y v 4 x 10*, for example. Note that if y F < 5 x 10 f ( n i / 

1 cm -3 )" 1/3 then the photon energies are below the pair pro- 
duction threshold (this condition is Lorentz-invariant, but is 
most easily evaluated using eq. [29] in the frame of the 
shocked fluid, in which the photon distribution is isotropic), 
and there can be no pair-production catastrophe. 

It has also been observed (Fishman 1993) that many GRBs, 
or subpulses within them, show a progressive spectral soften- 
ing with time. This is qualitatively explained using Figure l 
and equation (18). If the radiation field is isotropic m the frame 
of-fluid 2 (as will be the case if the particle distribution and 
magnetic field directions are isotropic) and has a characteristic 
photon energy, then in the local observer’s frame the spectral 
hardness above this characteristic spectral peak will be a 
decreasing function of 0, because of the angular dependence of 
the Doppler shift. Higher frequency photons are preferentially 
observed from smaller values of 6 , which arrive earlier in the 
burst or subpulse, while lower frequency photons are observed 
over a wider range of 6 and hence over a longer time. A quanti- 
tative prediction for the spectral evolution with time could be 
made by numerical integration of the synchrotron emission 
function, but would depend on (uncertain) assumptions made 
regarding the energy and angular distributions of the radiating 
electrons. Spectral softening within subpulses is also produced 
by the more rapid radiation of the most energetic electrons. 
The decline of y F as more interstellar matter is swept up leads 
to a spectrai softening throughout a GRB, from one subpulse 

When shock-accelerated relativistic electrons enter the inter- 
stellar medium they may coherently radiate synchrotron radi- 
ation in analogy to the coherent electromagnetic pulse (EMP) 
produced by terrestrial high-altitude nuclear explosions 
(Karzas & Latter 1962, 1965). Although radiation rates by 
individual electrons in interstellar fields are low, the coherent 
emission of many electrons may be significant. The observe 
synchrotron frequency of this natural EMP may extend up to 
the visible or ultraviolet bands, but the intensity may be 
reduced at higher frequencies at which the electrons do not 
radiate coherently. The cutoff - frequency depends on the elec- 
tron density as well as on their Lorentz factor, and is difficult 
to estimate. The problem is particularly complex because the 
debris energy density far exceeds that of the ambient field, 
which is strongly affected by the entry of the Par- 
ticles (even aside from the plasma processes in the . 

shock), ih contrast to the terrestrial case. At i^SaUctic 
an initial impulsive pulse is strongly affected by intergalactic 



f 


TT" 


plasma dispersion, with an electron column density 
— 10 23 (0_/0.1) electrons cm -2 , where is the fraction of the 
cosmological closure density present in ionized intergalactic 
gas. The resulting pulse may be very long (10 9 Hz radiation is 
delayed by a few minutes for flL = 0.1), and therefore difficult 
to detect. If it were observed it would have a characteristic 
dependence of frequency on time [v oc(t — *o) 1/2 > a cosmic 
whistler”, although in a terrestrial whistler v oc (f to) » Stix 
1962] which would directly determine the mean intergalactic 
electron density. Unfortunately, the limiting amplitude of the 
radiated EMP field is comparable to the ambient interstellar 
field in the source region (the EMP field screens the electron 
current from the ambient field, and limits^ the mdiated 
intensity), leading to an upper bound ~£ 2 SM cr 2 /2 ~ 10 ergs 
s” 1 on the radiated power which is very small. 
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3. POPULATION G: GALACTIC GAMMA-RAY BURSTERS 

GRBs which show spectral features, typically around a few 
tens of keV and at 400 keV, have long been identified with 
Galactic magnetic neutron stars and are inexplicable at cosmo- 
logical distances. Their distances cannot be determined from 
available data, and could be less than 100 pc, —100 kpc, or 
anything in between. The familiar arguments concerning the 
mechanisms of GRBs at Galactic distances center on two 
issues : the source of energy, and the physical conditions in the 
emitting region. The problems are harder, the greater the 
assumed distances. The observation of the 1979 March 5 event 
at a likely distance of 55 kpc (Cline 1980) forces the consider- 
ation of distances of that order, and of correspondingly high 
luminosities, even though it is unclear whether it (a SGR) was a 
member of the population G of GRBs, or represented a distinct 
third class of objects. 

The central problem of distant and luminous gamma-ray 
sources is gamma-gamma pair production (Cavallo & Rees 
1978; Schmidt 1978; Katz 1982; Epstein 1985; Carrigan & 

Katz 1992). This process does not permit the escape of a large 
luminosity of MeV gamma-rays from a small region unless 
they are collimated, and thus excludes many models of GRB or 
SGR at Galactic halo or cosmological distances. It is well 
known that this problem is avoided in a collimated relativistic 
outflow of radiating matter, a consideration which led to the 
popularity of fireball models, in which an opaque cloud of 
radiation and pair gas adiabatically expands and cools until its 
particles’ velocity vectors are collimated outward. However, 
fireballs are conspicuously incapable of producing low redshift p 
(400 keV) annihilation lines, line J^tures at tens of keV, or the T 
observed long and complex time structure. The interaction of 
fireballs with their environment may solve the temporal 
problem, as discussed in § 2, but offers no hope of solving the 
spectral problem. The case for magnetic neutron stars for 
GRBs or SGRs with spectral lines remains as strong as the 
data. 

It is usually assumed that the radiating region of a GRB or 
SGR in a neutron star model is dominated by pair plasma, 
with n + P where n + and n, are the positron and ion den- 
sities, respectively. This assumption is made, in analogy to 
fireballs, even for nonfireball models, perhaps because the 
threatened gamma-gamma pair production catastrophe seems 
a likely source of dense pair plasma, and because the observed 
annihilation line requires a source of positrons. However, the 
assumption of pair dominance may not be justified in non- 
fireball models, such as are required to explain population G 
GRBs. If sufficient gamma-ray collimation is present to avoid a 
gamma-gamma pair production catastrophe, then the pro- 
duction of pairs may be negligibly small When the observation 
of annihilation radiation provides empirical evidence for the 
production of some positrons, it should be remembered that an 
observably narrow annihilation line requires temperatures less 
than 50 keV, and may be produced by a comparatively small 
number of positrons precipitated onto the cool neutron star 
surface; a hot pair plasma does not produce a recognizable 
annihilation line. 

If a pair plasma is not an expanding fireball, it must be 
trapped on magnetic field lines (Carrigan & Katz 1992). Gravi- 
tation is unimportant for pairs, so they fill a magnetosphere 
(presumably of a magnetic neutron star). However, they 
quickly (in a free-flight time) precipitate onto the stellar 
surface, where they annihilate, because they more rapidly 
radiate their transverse momentum^ the cyclotron process; 
even if the radiation density is4*msufflcient to maintain most 
leptons in excited Landau (magnetic) states (a condition satis- 
fied under only the most extreme conditions), their interaction 
with the radiation field destroys their transverse adiabatic 
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invariant. In this magnetosphere-filling geometry the emergent 
radiation, by whatever process, is not collimated, and gamma- 
gamma pair production imposes its usual limits on the emer- 
gent flux of Mev gamma-rays. 

It may be more satisfactory to consider an electron-ion 
plasma with only a small admixture of positrons (sufficient to 
produce the observed annihilationjine), so that n+ <$ n ( . An 
electric field may accelerate jwitnthe electrons, which radiate 
by bremsstrahlung or by the cyclotron process after elastic 
scattering on the ions raises them to excited Landau states. ^ 
Because most of the leptons are negative, they form a broadly 
collimated beam; if they are relativistic the resulting radiation 
is similarly collimated and there is no gamma-gamma pair 
production catastrophe or limit (other than a Plarjc function at 
an effective temperature characterizing the elecffon distribu- 
tion function) on the emergent intensity. 

In contrast, an electric field acting on a pair gas heats it but 
imparts no net momentum to the leptons; two counter- 
streaming beams of gamma-rays readily produce pairs, rapidly 
achieving equilibrium with them and limiting the emergent 
intensity. A minority admixture of positrons in an electron-ion 
plasma produces only a proportionately small counter current 
of gamma-rays to those produced by the electrons. This 
countercurrent removes an equal current of electron-produced 
gamma-rays by pair production, but the remaining electron- 
produced gamma-rays form a collimated beam and escape 
comparatively freely, suffering little or no (depending on the 
degree of collimation) gamma-gamma pair production. This is 
a consequence of the net momentum imparted to the lepton- 
photon system by the electric field in analogy to the momen- 
tum imparted to a sector of a fireball by adiabatic expansion. 

It is possible to make simple rough estimates of the param- 
eters of an electrically heated ion-electron sheet plasma, which 
might be the source region of a GRB of populationG, follow- 
ing Katz (1993). Consider a sheet of thickness L,-becomposed 
of positive ions of charge Z and density (n e = Zn,), .have 
transverse optical depth t and temperature T, and radiatih^ 
power per unit area P. Define the dimensionless power per unit 
area p = Ph* /(m* c% dimensionless thickness / = 
Lm e c 2 /e 2 , and temperature t s k B T/{m e c 2 ), where these quan- 
tities have been scaled to values characteristic of a relativistic 
electron (or pair) gas. The characteristic radiant intensity 
m*c 6 /h* = 4.3 x 10 35 ergs/(cm 2 s) and length e 2 ^/ 2 ) = 2.8 
x 10" 13 cm (the classical electron radius). The optical depth is 


r ~n e a 0 L, (31) 

where the characteristic cross section <r 0 = e 4 /(m 2 c 4 ). At non- 
relativistic energies the appropriate cross section is the 
Thomson cross section 87rcr 0 /3, but at semi-relativistic energies 
of interest a 0 may be a fair approximation. The observation of 
a nonthermal spectrum implies that t cannot much exceed 
unity, but t ~ 1 and t 1 are each possible. 

The large field of a magnetized neutron star has a number of 
effects. It enters the argument of the effective Coulomb 
logarithm in collisional processes, typically reducing it to 
In A % In [fc B Tm e c/(heB) ] (Katz 1982). Both bulk motion and 
current flow are restricted to be parallel to the field lines, jus- 
tifying the assumption of thin sheet geometry and making the 
field distribution nearly force-free (J x B = 0). Perhaps most 
important, it means that any electron energy resulting from 
motion perpendicular to the field is immediately radiated. 
Even in conditions characteristic of the 1979 March 5 event at 
55 kpc the radiation density is far below Planckian (for f ~ 1), 
so that electrons may be assumed to be in their ground mag- 
netic state until collisionally excited, and then to radiate as if in 
vacuum. The radiation rate per unit area may therefore be 
estimated using standard expressions for elastic scattering: 


P ~ n e nik B T 
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This expression may be inverted, using the definitions of p, i, r, 
and a = e 2 /(hc\ to give 


pt 1/2 V /2 m 3 c 3 


\la*Zln\J ft 3 ’ 

the characteristic density m 3 c 3 /fr 3 = 1.8 x 10 31 cm 3 . An 
alternative expression for n e is obtained from equation (31): 


n '~/a 3 h l ' 

Equating the expressions (33) and (34) yields a result for the 
thickness: 

l ~^nrr- ( 35 > 

pt 112 a 3 

Very roughly, p ~ 10' 4 for the 1979 March 5 event, if in the 
LMC, and p ~ 10” 7 for a typical observed GRB if at 100 kpc 
distance, correspnding to L ~ 1.4t 2 cm and L ~ 1.4 x 10 t 
cm, respectively, where Z — 26, t = 0 . 1 , and A = 10 were 
taken. The densities are correspondingly high. 

If the radiating sheet is driven by magnetic reconnection, as 
is plausible, then the radiated power may be related to the 
electrical work done : 


where E is the electric field (properly, its component parallel to 
B) and a nonrelativistic expression (Spitzer 1962) is used for the 
(cgs) electrical conductivity : 




2 (fc„r ) 3/2 

m, l/2 e 2 fZlnA 
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where the parameter <J ^ 1 is a correction factor which allows 
for the possibility of anomalous (plasma instability) resistivity 
when current densities and electron drift velocities are large 
and for the decrease in conductivity when electron velocities 
approach c. The characteristic conductivity m,c 3 /e 2 =11 
x 10 23 s' l . Equations (36) and (37) may be combined with the 
definitions of the varous dimensionless parameters to give the 
electric field, current density, and electron drift velocity : 

( 38 ) 


j cf ‘ <2 Ah T\ il2 1 (40\ 

( ) 

The characteristic current density m 3 c 7 /e 5 = 6.5 x 10 3B esu/ 
(cm 2 s). The drift velocity is thus comparable to the electron 
thermal velocity unless £ P 1 ; ion-acoustic instability is likely 
unless f > mJ(Zm t ). 

GRB mechanisms based on magnetic reconnection, as in 
solar flares, have been qualitatively discussed for many years 
(Ruderman 1975). They may explain GRB energetics and 
phenomenology (Katz 1982). Magnetic reconnection by sheet 
currents provides a natural explanation of the electrically 
heated sheets discussed in this section. If this model is assumed, 
then Maxwell’s equations provide an additional relation 
among j , L, and the magnetic field B, permitting further con- 
straints to be placed on the parameters. Because j is parallel to 
B, the direction of B rotates across a sheet current without 
changing its magnitude B. If the total angle of rotation across a 
uniform current sheet is it, then 

Be 

L ~ 4;' 


(41) 
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Defining the usual characteristic magnetic field B c - 
m 1 c 3 /(eh) = 4.4 x 10 13 gauss and the dimensionless parameter 
b s B/{4B C ), equation (41) may be rewritten 


Tb£ 1/2 ZlnA 
pta 2 


Equating tty^expression to equation (35) yields 

bat 112 




(42) 


(43) 


For plausible t ~ 1 and b ~ 0.02, t ~ lO'^ 1 ' 2 ; either the 
emission region is very optically thin or the ^ivrt^is domi- 
nated by plasma wave scattering, and is far in excess oti 
independent particle value. Either or both of these possibilities 

1S The P electnca^eld (eq. [38]) may be evaluat^. ujng equa- 
tion (34) for n e and equation (35) for l, and ^fining the charac- 
teristic electric field E c = m) c 3 /{eh) m B c = 4.4 x 10 cgs . 
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If the energy release is driven by magnetic reconnection it is 
proper to use equation (43) for r, yielding 

(45) 


r P<X F 

E ~J E '- 


At Galactic halo distances the brightest GRBs may have E 
sufficient to produce vacuum breakdown into : a i paw ' 8“ (Srro th 
& Epstein 1993), but in less intense or closer GRBs this will not 
occur and the resistively heated ion-electron p asma d,scus«d 
here may be sufficient. On a microscopic level, of ^ our ^’ * 
power is determined by the magnetic field strength and con- 
figuration, and by the mechanisms of plasma resistivity which 

drive reconnection. . rPR 

The observed nonthermal gamma-ray spectrum of GRBs 
requires the presence of a nonthermal distribution of 
“runaway” electrons, consistent with the large v ir and collec- 
ti^ processes discus^d above. It may be relevant that values 
of r ~ 1, consistent with the radiation of the bulk of the power 
of GRBs at ~Mev energies, correspond to a maximum in the 
conductivity and therefore, under conditions of a fixed poten- 
tial drop, to a maximum in the power dissipation. . , 

If z P 1, cr t | increases approximately linearly with w+ in 
range Zn <n + < Z 2 n, because of the increasing density of 
charge carriers without a corresponding decrease m their sca - 
tering length.TKis is in contrast to the usual near-independent 
of density of <r el . Pair production thus may provide a natural 
thermostat at t ~ 1 reducing the power dissipated in regimes 
it the pair production threshold by increasing the conductivity 
under coitions of constant current. It is evident, of court 
that observable cyclotron and annihilation 1 mes require muc 
lower values of t, and are plausibly produced by energy and 
positrons precipitated on the dense surface layers of the 
neutron star, cooled by blackbody radiation. 

4. DISCUSSION 

The fundamental problem of 

apparent inconsistency between their spatial distnbut'on 
which points strongly toward a cosmological ongm, and n “* 
spectral features observed in some GRBs, apparently incons 
tKth such an origin. -In this paper, I have tned to reconcile 
these apparently contradictory data by assuming two disjoint 
populations of GRBs. Because the argument for a cosmo- 
EScal population C is statistical, while the argument for a 
Galactic population G is based on the observation of s P«ctr 
lines from only a minority of GRBs (perhaps from only a 
minority of that population), it is difficult to assign an ind.vid- 
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ual GRB to either population, unless it shows spectral lines or 
is identified with another astronomical object (the SGR of 1979 
March 5 is the only good extant example of such an 

identification). ... . 

Fortunately, the models discussed in this paper predict 
another potential distinguishing characteristic. GRBs in popu- 
lation C produce their radiation by the synchrotron process of 
relativistic electrons, if the magnetic field is ordered the radi- 
ation will be linearly polarized. GRBs in population G produce 
their radiation by the cyclotron process of semi-relativistic 
electrons. Coulomb scattered into excited Landau states. This 
radiation is elliptically polarized, with a substantial circular 
component. In contrast, radiation produced by annihilation m 
a pair gas, such as the initial burst from a fireball or the 
trapped pair gas discussed by Carrigan & Katz (1992), is unpo- 
larized. 

The predicted characteristic frequency of gamma-ray emis- 
sion in fireball debris impact models (eq. [30]) is fairly sensitive 
to the initial fireball Lorentz factor y F . The value of y F depends 
on physical conditions within the fireball (Shemi & Piran 
1990), and a wide range of y F would be expected, with a corre- 
sponding range in spectra. In particular, smaller values of y F 
would lead to X-ray, ultraviolet, or visible bursts without 
gamma-ray emission; these should be searched for. Bursts of 
lower frequency radiation should have longer durations and 
smoother time histories (eq. [28]), and perhaps also lower effi- 
ciencies, as accelerated electrons undergo adiabatic expansion 
(coupling their kinetic energy to ion motion) before they 

Relativistic synchrotron models, such as those discussed 
here, generally involve a power-law distribution of radiating 
electrons up to a cutoff Lorentz factor y 2< (eq. [14]), and there- 
fore a power-law synchrotron spectrum up to a cutoff fre- 
quency v 2 (eq. [29]). The power-law indices are uncertain, but 
using the usual estimate for the spectral index p = (r? + 2)/ 

(q _ 1) of relativistic particles accelerated by a shock of com- 
pession ratio rj and the relation s = (p — l)/2 for the 
synchrotron spectral index yields s = 3/[2(^ — V )] = 
3(y fl - l)/ 4 = i for a strong shock in a relativistic fluid for 
which ri = l and the adiabatic exponent y a = 4/3. Because 
p = 3/2 most of the particle energy resides in the most ener- 
getic particles, in contrast to the usual case with p > 2. The 
radiation spectrum F v oc implies a visible power - 10' of 
that in soft gamma-rays; a bright GRB of 10' ergs/(cm s) 
extrapolates toF,-5x 10' 23 (v/10 9 Hz)' 1/4 ergs/(cm 2 s Hz). 
However self-absorption limits the brightness temperature at 
radio frequencies. When the debris shell has a radius r, the 
effective radiating area is ~ r\l y F bqpyuse of the collimation, 
and the brightness temperature T b < y F rn P c Ab* Th e resulting 
power spectral density at a distance D is < 2tcv tn p r t /D ^ 
x 10' 29 (v/10 10 Hz) 2 (rJ2 x 10 15 cm) 2 ergs/(cm 2 s Hz). As the 
debris shell expands this may become significant. Lower 
energy electrons radiate more slowly in a given field, with r < * 
v “ */* (in any frame), so that an observed t r of 1 msec at 300 kev 
corrsponds to a minimum pulse length of -0.3 s in visible light 
and of - 10 2 s at 1 GHz. Radiation from expanded shells with 
r, > r,, where y F is too small to produce gamma-rays and B is 
small, will have much longer pulse lengths. 

Relativistic electrons of lesser energy will also be produced 
after the debris have been slowed by interaction with sur- 
rounding gas (and gamma-ray emitting electrons are no longer 
accelerated), so that visible and radio-frequency pulses may 
follow the high energy pulse by seconds, minutes, or even 
hours, and be characterized by smaller y F and correspondingly 
longer durations. 

During the lengthy process of review of this paper several 
groups (KouveliotO-^t al. 1993; Lamb ly Graziam,& Smith 
1993- and Mao, Narayan, & Piran 199?) reported that the 
durations of GRB are bimodally distributed. This may be 
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explained, using equation (28), as a consequence of the well- 
known bimodal (or multi-modal) distribution of interstellar 
gas densities in ionized and neutral regions; the observed dura- 
tion of a GRB is determined by y F (which may be distributed 
over a wide range) and the gas density in a small region (see eq. 
[2]) surrounding the fireball. 

The natural next step in the development of population C 
GRB models described here would be to perform numerical 
one-dimensional special-relativistic hydrodynamic calcu- 
lations of relativistic debris clouds propagating into model dis- 
tributions of clumpy interstellar media. The predicted pulse 
shapes and spectra may then be calculated using equation (24) 
convolved with the radiation timescale of electrons as a func- 
tion of their energy. 

I thank P. C. Joss and I. A. Smith for discussions and NASA 
NAGW-2918 for support. 




APPENDIX 


As relativistic debris sweeps up surrounding diffuse matter it forms a relativistic blast wave. The debris energy degrades. Its 
radiation shifts to lower frequencies and increases in duration, and may be observable at frequencies from X-rays down to radio. 
This appendix collects a few simple results concerning relativistic blast waves. 

First, assume that a negligible fraction of the blast wave’s energy is radiated, as in the Taylor-Sedov nonrelativistic blast wave. As 
long as the blast wave remains relativistic its radius r % ct, where r and t are its radius and age as measured by a local observer. The 
dependence of other quantities on r and t is more interesting. Conservation of energy gives 


yA r ) * 


3 Y 
An r z pc 2 


(46) 


Define a characteristic radius 



(47) 


where y 0 s y^O) is the original debris Lorentz factor, so that 

fcW-7^". < 48 > 

The relativistic particles, after being thermalized in a shock, are roughly isotropic in the frame of that shock and have an angular 
dispersion 9 0 ~ }y 1/2 (eq. [ 1 7]) as measured in the local observer’s frame. As a result, they fill a shell of thickness 


Ar 



(49)' 


where the chief contribution to Ar is their angular dispersion (their dispersion in speed contributes only a thickness 0[r/y}]). The 
energy density measured in the local observers frame is then 


Y 3 Y 1 

<?(f) ~ 4nr 2 Ar ~ 32* 2 r 6 pc 2 ' 


(50) 


while the proper energy density is less by two factors of y F /2 : 


e(r) 


Y 
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(using eq. [46], this is equivalent to eq. [10], up to factors 0[1]). Then, as in equation (26), the magnetic field in the comoving frame 
is 



(52) 


In this frame the electric fields are small and the usual synchrotron radiation expressions may be used (in the local observer s frame 
there are large electric fields which reduce the electrons* magnetic acceleration and permit them to remain collimated despite the 

presence of transverse magnetic fields). , 

Using the previous results for electron energies, the characteristic synchrotron radiation frequency, measured in the comoving 
frame, is 


v 2 (r) 9 2 


hiA. ( 53 ) 
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where v 2 (r 0 ) is given in equation (29). In the local observer’s frame this frequency is 

v ob,M ~ v ob ,(r 0 )^j , 

where v obs (r 0 ) is given in equation (30). The electron energy-loss time, measured in the comoving frame, is 

'lo» ~ tr2(r 0 )(^ , 

where t rl is given in equation (27). In the local observer’s frame this time is 

b,M ~ . 

where t obt (r 0 ) is given in equation (28). 

The radiated power, measured in the comoving frame, is 
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while the radiated power measured in the local observer’s frame is 

2n Ze*Yy? 
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The efficiency of radiation over a range ~ r around radius r is 

P 0 Jr) r 2n {e*y% 2 


Q(r) ~ 
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the factor y\ l2 (r) results from the Doppler shortening of the observed duration of the radiation. 
Defining a characteristic time 
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where f 0 = r 0 /( r >’o ,2 )> yields the scaling laws : 
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PcJtc ) ~ P -(t 0 )(^) 6/5 , 
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Numerical evaluation shows that e(r 0 ) |> 1 for plausible parameters. Such a high nominal efficiency of radiation justifies the 
assumption of efficient gamma- ray production, but is inconsistent with the assumption of energy conservation made to obtain the 
blast wave results of equations (46)-(69). The opposite limit, in which the collision between fireball debris and interstellar matter is 
completely inelastic (analogous to the late-time “ snowplow ” limit of nonrelativistic blast wave theory) is easily treated. 

In a narrow angular segment dQ a debris proper mass MdCl/(4n) = Y dfi/(47ry 0 c 2 ) sweeps up an interstellar proper mass pr 2 dr dQ 
in traveling a thickness dr . In the local observer’s frame the interstellar matter is at rest, and the momentum of the debris is 
Mc(y 2 - l) ij2 dCl/(4n). After the collision the total proper mass is [M/(4 tt) + pr 2 dr]dCl. The Lorentz factor y + dy of the new 
center-of-momentum frame is found from the expression p. if a relati vistically moving object of the total proper mass : 

hr 
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dy 2 - 1) 1/2 = 
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+ pr 2 dr dnjc[(y + dy) 2 — 1] 1/2 . 


(70) 


If the radiation is emitted isotropically in the center-of-momentum 
Zlementary calculation yields the differential equation 


frame its emission does not change the velocity of this frame. 


dy _ 4nr 2 p y 2 — 1 
dr My 


(71) 


tegration yields the result 


y 2 - 1 = (y 2 o - 1) exp [-87ur 3 p/(3M)] 

= (Vo ~ 1) exp (—2 r 3 /rl) . (72) 

The relativistic debris are rapidly degraded by radiation at r a r 0 . Once y decreases significantly the efficiency of radiation 
dines and the non-radiating blast wave results (46)— (69) become applicable, with Y, y 0 , and r 0 assuming the values appropriate to 
5 debris at the transition between the efficiently radiating and weakly radiating regimes. In practice, while gamma-ray emission 
nay occur in the strongly radiating regime, emission of lower energy photons will be predominantly from the weakly radiating 
gimes. 

The exponential dependence of y on r 3 (eq. [72]) implies that the transition between the two regimes occurs for r ~ O(r 0 ), even 
ough the degradation of Y and y may be large. This degradation may be estimated by setting e(r) = 1 and adopting r a r 0 . The 
rviving fireball energy Y' and Lorentz factor y' may be approximated 

r ~ re 0 “ 2/3 , 

y'o ~ 7o e o 2/3 . (73) 

here e 0 is the nominal radiative efficiency obtained from equation (59) using the original fireball Y andy 0 . The surviving Y' and y' 0 
ay then be used to describe the subsequent evolution and radiation of the blast wave. 
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Note added in proof: 

The analysis of the flow shown in Figure 1 may readily be extended to the case n \ ^ ?i 4 . 
Then e-i = e 2 = 27F(e ] e 4 ) 1 / 2 = 27/rm a c 2 (ni?r 4 ) 1 / 2 . The Lorentz factor describing the 
speed of the shocked matter in the local observer’s frame 712 = (7F/2) 1/,2 (n 4 /ni) 1 / 4 . The 
internal energy of the shocked interstellar matter is described by the Lorentz factor (eq. 
[ 13 ]) 72 = (7F/2) 1//2 (n4/ni) 1/f4 . The shock compression n 2 /nj = 4712. The equipartition 
magnetic field B 2 oc t^ 2 oc (nin 4 ) 1/ ' 4 , so that the co- moving synchrotron frequency v 2 cc 
72B2 c < (n 4 /ni) 3 / 4 , and the observed frequency u 0 b s oc n 4 /n\. Unshocked relativistic shells 
axe very thin (A rjr ~ 7p 2 ) so that n 4 n 1 is expected; this substantially reduces the 7 f 
required to obtain the observed high energy gamma-rays. 

The passage of the shock S2 through the debris shell leads to a reflected rarefaction 
which reduces both n 4 and 7 f- In the local observer’s frame us2/c ~ 1 — (2/7/r)(?rj /n 4 ) J ' 2 
(eq. [ 7 ]). If n 4 /ri] 0 ( 7 ^) then 1 - VS2/C ^ 0(7 r 2 ) and S2 will not pass entirely through 
the shell until Si has passed through a thickness 0 (r) of interstellar cloud. In contrast, 
for TI4 ~ Tti S 2 traverses the entire debris shell and is reflected from its back surface as 
a rarefaction after Si has traveled only 0 (r/ jf)- Qualitatively similar remarks apply to 
debris shells with continuous distributions of 7/^. 

T. Piran has independently obtained similar results. 



